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CO ■ Abstract 

>. 

C*") ' An attempt is made to go beyond the semi-classical approximation for gravity in the 

Born-Oppenheimer decomposition of the wave-function in minisuperspace. New terms are 
included which correspond to quantum gravitational fluctuations on the background metric. 
They induce a back-reaction on the semi-classical background and can lead to the avoidance 
QO \ of the singularities the classical theory predicts in cosmology and in the gravitational collapse 

0^ ■ of compact objects. 
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1 Introduction 



The canonical quantization of highly symmetrical general relativistic systems carried out in 
suitably chosen variables leads to the dynamics being determined by the (super)Hamiltonian 
constraint [jl], |2| of the Arnowitt, Deser and Misner (ADM) construction || in the space of 
functions of time called minisuperspace [|l], |[]. Such an approach is particularly useful to inves- 
tigate self-gravitating quantized matter with gravity in the semi-classical regime. One performs 
a Born-Oppenheimer (BO) decomposition of the wave-function satisfying the Wheeler-DeWitt 
(WDW) equation into two parts ||. The first one represents a collective degree of freedom 
associated with gravity (slow component) and, in the semi-classical approximation, leads to an 
Hamilton-Jacobi (HJ) equation for the gravitational degree of freedom; the second part describes 
microscopic matter (fast component) and satisfies a Schrodinger equation in the time defined by 
semi-classical gravity. 

Alternative approaches have been attempted, such as the one in Ref. ||, where, however, 
relative phases of matter and gravity were incorrectly identified fH , or which involve an expansion 
in the Planck mass (see, e.g., Ref. 0, ||). The latter expansion is potentially dangerous, since 
it has been shown that it can lead to violation of unitarity within the framework of canonical 
quantization ||] and to incorrect identification of the background as an empty solution of Einstein 



equations [||, 10]. In the BO approach M the collective degree of freedom evolves slowly because 
it is associated with the total mass of the system which is (many) times the mass of each 
constituting matter quantum (regardless of the latter being smaller than the Planck mass). 
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For these reasons we shall appeal to the approach introduced in Ref. || as the best suited for 
the purpose of analyzing the semi-classical limit and shall not attempt at expanding the total 
wave-function in the Planck mass. 

The BO approach was subsequently applied to two physical models of general interest: the 
gravitational collapse of a sphere of homogeneous dust in empty space (TJ| and spatially homo- 
geneous Universes |9|, ^] (see also Ref. [13] for collapsing shells). For the former system the 



novel effect of non-adiabatic production of matter has been studied with the analytical method 
of the (adiabatic) invariants for time dependent Hamiltonians |l4], ||] in Ref. [15|. The same 



technique, supplemented by numerical simulations, has shown the possibility of having an infla- 
tionary phase in the primordial Universe which is driven by purely quantum fluctuations of the 
inflaton and has finite duration [jilf . In both cases there are one degree of freedom for gravity, 
R (related to the external radius of the sphere or the scale factor of the Universe), and one 
degree of freedom for matter, (ft (homogeneous scalar field). The phase space is then the usual 
Friedmann-Robertson- Walker (FRW) minisuperspace of the space-time metric 



ds 2 = R 
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(e = 0, ±1 respectively for flat, spherical and hyperbolic space and p > for e = 0, —1; < p < 1 
for e = 1) and homogeneous scalar matter. The coordinates (rj, p, 9, ip) define a comoving 
reference frame and for the sphere of dust p < p s , where p s is the (constant) comoving radius 
of the sphere. 

One must be careful in modeling dust with a scalar field, since the latter indeed describes a 
perfect fluid with pressure equal to the Lagrangian density, p ~ \ ((ft eft 2 ) WA - If the scalar 

field has mass = h/i^, then p oscillates with frequency ~ 2/1^, e.g., for ~ 10~ 27 kg (the 
proton mass) this means a period T ~ 1CP 23 s. It is thus reasonable to approximate the actual 
pressure with its time average over one period (that is, set p = 0) provided the radius R does 
not change appreciably on the time scale T (quantum adiabatic approximation for the state of 
the scalar field). Moreover, this adiabatic approximation becomes exact in the classical limit 
for (ft, as can be seen by taking ft — > with held fixed (£^ — > and T ~ vanishes), and a 
mode of the homogeneous massive scalar field can be identified with dust. In fact in Ref. [^TJ 
it was verified that in this approximation one recovers the classical Oppenheimer-Snyder (OS) 
model Q. 

A major restriction in Refs. H [9|, 0, [l^, 13, [H| 16] is that quantum fluctuations of the 



gravitational degree of freedom were suppressed a priori and R was approximated by a classical 
trajectory i? c (r/ c ) (rj c being the conformal time associated with that trajectory) which, in turn, 
was determined solely by the matter content. The aim of the present notes is to allow the variable 
R to have quantum fluctuations around the classical trajectory and modify the expressions of 
the general formalism || accordingly. Of course, it would be much more interesting to allow for 
inhomogeneous fluctuations, but this would inevitably render the system intractable analitically 
and is left for future developments. As a by-product, we will see that one can have a significant 
back-reaction of the gravitational quantum dynamics on the semi-classical trajectory. This 
affects the singularity classical General Relativity generically predicts in cosmology and as the 



final state of a collapsing body (see IE] and Refs. therein). 

The possibility of avoiding space-time singularities in a quantum theory has been studied for 
a long time and the literature on this topic is wide. Here we only refer to two approaches: 
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1. in quantum field theory in curved space-time (see, e.g., Ref. p0| ) gravity is described by a 
classical background on which quantum matter fields propagate. In Ref. it was found 
that there are states of matter for which the Universe admits a minimum non-zero scale 
factor, provided the number of particles is not conserved; 

2. it was suggested that canonical quantization of the gravitational degrees of freedom could 
bypass the cosmological singularity (2^]. In Ref. [^3[ the constraints were implemented 
before quantizing and one ended up with quantized gravitational degrees of freedom only. 
In this case no significant change in the classical behaviour was found. 

It is a trivial observation that in a quantum theory a point-like singularity is meaningless since 
it would violate Heisenberg's principle. What we shall show in the proposed approach is that 
one expects the singularity is avoided under a broad assignment of initial conditions. In fact the 
semi-classical approximation breaks down before the point-like singularity is reached (but within 
the adiabatic approximation for the gravitational degree of freedom) and the very concept of a 
trajectory loses its meaning at a value of R which can be appreciably big (in a sense that will 
be specified later). 

The plan of the paper is as follows. In the next Section quantum gravitational fluctuations 
are treated in the standard BO formalism for the FRW minisuperspace and it is shown that 
their energy cannot always be neglected with respect to the energy of matter. In Section ||| 
the energy of such fluctuations is incorporated in a modified semi-classical HJ equation which 
is solved under certain approximations. Such approximations are then analyzed to determine 
the range of validity of the solutions. In Section || some conclusions are drawn for cosmological 
models and for the collapse of homogeneous spheres of dust. Finally in Section ||| the results are 
summarized and commented, we shall use units in which c = 1, k = 8ir Gn, £ p = VhK is the 
Planck length. 



2 Quantum gravitational fluctuations in the BO approach 



Let us start directly from the WDW equation in the minisuperspace of the two variables R 
and 4> (f° r a derivation from first principles see [|J) with a convenient operator ordering in the 
gravitational kinetic term [11]: 



Hq + Hm 
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The wave- function ^ can be expressed in the factorized form ^{R 



V(R, 







(2.1) 



Rip(R)x(<l>,R) which, 



after multiplying on the LHS of Eq. (2.1) by x* an d integrating over the matter degrees of 
freedom, leads to the equation for the gravitational part B 
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The scalar product (x\x) — I d4>x*(4>,R)x(4>>R) an d 



7/} = e -iJ R A(R')dR'j 



X 



,+t f R A(R')dR> ~ 



X 



(2.2) 



(2.3) 



3 



with A = — i ( x | X ) 1 ( X I 9r \ x ) = ~ * ( 9r ) ■ If we now multiply Eq. ( |2.2|) by x and subtract 
it from Eq. fl2.1|) we obtain the equation for the matter function x 



H M -(H M ) x + Kh 




k% 2 



d 2 

dR 2 



d 2 

dR 2 



X 



(2.4) 



The Eqs. ( |2.2j ) and fl2.4j ), as well as the WDW equation, are exact, in the sense that no approx- 
imation has been assumed yet for the wave-functions x an d tp, and contain no time variable. 

A way one can introduce the time is by taking the semi-classical limit for gravity [^4|, B5l |5| . 
In order to do so, first one needs to neglect the RHS's of Eqs. ( |2.2| ) and ( p.4[ ) which are related 
to quantum transitions among different semi-classical trajectories [[§. As usual || [J, we shall 
check the consistency of all approximations once the solutions to the semi-classical equations 
have been obtained (see Section |3.2|) . In fact, it is not necessary (nor possible, in general) to 
prove that the RHS's are small from the onset, but it is sufficient to show a posteriori that they 
are negligible for the cases considered. Then one writes a semi-classical (WKB) approximation 
for the wave function ip 



i 



PcdR c 



where 



Pr 



1 dR, 



k drj c 



2kR c (H m ) -eR 2 



(2.5) 



(2.6) 



is the canonical momentum conjugated to R in the classically allowed region ( Hm ) > eR 2 /2k 
and the integral in the exponent is computed along the (so far unspecified) semi-classical tra- 
jectory R = R c {ric) with momentum P = P c (r] c ). Moreover, the derivatives with respect to the 
conformal time r\ c are defined according to Eq. (|2.6| ) as 



9 , p 9 8 

« = —KW r P— — = —KP r — — 

dr] c Vc dR dR 



Rc 



(2.7) 



where the last step follows from ip c having support only for R ~ R c , rj ~ rj c . 

Upon substituting = ip c into Eq. (|2.2| ), the gravitational equation finally reduces to the 
semi-classical HJ equation 



1 fdR\ 2 € 



2k \dr) ) 



2k 



R 2 + R{H M , 



1 fdRc 



2k \ dr] c 



R 2 +R C {H M ] 

2k 







(2. 



which can now be used to determine R c explicitly once ( Hm ) is given. 

It is important to note that the semi-classical regime is not defined simply as the limit fi — > 0, 
but rather by a specific choice of the wave-function ip c . For instance, with P c given by Eq. ([2.6]), 
( Hm ) = and a constant radial number density of scalar quanta > (in practice this 

is the statement of the quantum adiabatic approximation), one has (see, e.g., Ref. p6[) 



R, 
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(2.9) 
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that is the usual FRW cosmological models for increasing rj c or the OS model of gravitational 
collapse for decreasing r\ c (the classical singularity occurs at r\ c = in both cases). 
Substituting ip = ip c in Eq. (|2.4|) gives the Schrodinger equation 



i h 



dr] c 



4 2 



Xs 



(2.10) 



for the rescaled matter function x = Xs exp{(i/ft) f£ R c dr]' c ( Hm )}• We note in passing that the 
difference between Xs and the original x amounts exactly to the phase factor between eigenvalues 
of an hermitian invariant for the Hamiltonian in Eq. ( j2,10|) 14] and exact solutions of Eq. ( |2.10| ) 
pi. The above Schrodinger equation together with the HJ equation ( |2.§| ) was the starting point 
for the results found in Refs. [11, 12, 13, 15, 16] and led to the conclusions briefly mentioned in 
the introduction. 

So far matter and gravity are determined by two equations of clearly different types. Suppose 
instead one defines 



*l> = Af, (2.H) 
where / = f(R) will encode quantum fluctuations around the trajectory R c swept by ip c . Then 



Eq. fl2.2|) , again neglecting the RHS's (see Section |3_^ for a detailed description of the difference 
with the previous case), becomes an equation for f ( ' = d/dR) 
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(2.12) 



Upon using the definition ( |2.7| ), the above takes the form of a time-dependent Schrodinger 
equation ( ' = d/drj), 
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dr]c 
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(2.13) 



In the RHS, due to the factor tp c , the quantity P is still evaluated at the classical momentum 
and r] = rj c . Therefore one can simplify Eq. ( 2.13p by making use of the HJ equation (2.8) and 
obtains 



i h 



df_ 

dr/c 



2 



Ak 2 P 4 2k 2 R3 J W t K P 2 dR 



d 2 f 

R C + dR 2 



(2.14) 



which is now the analogue of Eq. ( |2.10| ) for gravity. 

Since Eq. (|2.14| ) is involved, let us take the classical adiabatic limit for R c , to wit |i? c | <C 
R c \d n+1 R c / drfc +l \ <C |i? c | for any integer n > 0. This approximation is not to be confused 
with the previously mentioned quantum adiabatic limit on the state of 4> an d will be further 
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discussed in Section |3.2| . From the definitions (|2,6|) and (2.7) it then follows that one can neglect 
terms containing P c and P c . Hence, only the last term survives in the RHS above and one finds 



df nh 2 d 2 f 

i n 



dr] c 2 dR 2 



(2.15) 

Rc 



which resembles the non-relativistic equation for a free particle of "mass" 1/k and negative 
kinetic energy. 

Plane waves 

Eq. ( 2,15| ) admits solutions in the form of plane waves, 



/A = eXP {*f5 + *f} ' (2 - 16) 

where the A's are real numbers (A > for left movers and A < for right movers in R space). 
One observes that, due to the "wrong" sign mentioned above, the energy (conjugated to the 
proper time dr = R c drj c ) associated with each mode of wavelength A, 

is negative. Although disturbing at first sight, this is in agreement with gravity contributing 



negative amounts to the total (super) Hamiltonian 27]. 



Exponential waves 



For imaginary A = il, from Eq. ( 2.16 ) one obtains a new set of solutions given by 



whose energies are positive, 



Of course the amplitude of the above solutions increase with R c for I > and decrease with R c 
for I < 0, thus signalling an instability. In fact, these modes can be related to the tunneling of 
R across classically forbidden regions, since iP (In//)' is real if P is imaginary. For this reason, 
one cannot superpose solutions from the two sets {fx} and {/;} if the classical limit for R has 
to make sense |28| (due to the factor ip c either R is in a classically allowed region or it is not). 
In the following we will only consider real values for A. 
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Full solutions 



The full gravitational state corresponding to the modes f\ found above are given by 

VA^^/A = ^exp|i|j + ^| , (2.20) 

where the weight ip c ensures that / dRip^O ip\ = (O f\)(R c ) for every operator 0(R, d/dR). 
Then the general solution to Eq. ( 2.15| ) is a superposition of the form J2 c\ f\{R c ) and the total 



energy associated to quantum gravitational fluctuations is given by 



h£ 2 \r^\ 2 
c A 



where the c\ are normalization coefficients. 



It is clear from Eq. (2.21) that \Ef \ can be very large, depending on the modes fx which are 



included. In particular, \E\\ <C (Hm) only for 

f, t 2 

A 2 » \\ = J %— . (2.22) 

2Rc(H M ) 

An interesting observation is that A c is time-dependent (via R c ) and, in the quantum adiabatic 
approximation for Xs, the term (Hm ) = N^m^ is constant |ll[] and one has 

For an expanding universe in which R c increases in time without bounds, A c will eventually 
vanish after it had been as big as possible in the far past. On the other hand, for the case of 
a collapsing sphere of dust with monotonically decreasing R c , A c will diverge and, no matter 
how long are the wavelengths of the initial gravitational fluctuations, \Ef\ will overcome ( Hm ) 
before the sphere reaches the classical singularity R c = 0. 

We thus arrive at the following paradoxical conclusion. Our equations show that there are 
quantum gravitational fluctuations which can be generally associated to the classical solutions 
R c in Eq. ( |2.9[ ) . The energy of such fluctuations becomes inevitably larger then the matter energy 
at certain times but, since it does not appear in the HJ equation, the presence of gravitational 
fluctuations does not affect the semi-classical motion in any way. In the next Section we shall 
show that this paradox is due to an incorrect identification of the semi-classical limit and propose 
an approach to include the back-reaction of gravitational fluctuations. 



3 Improved BO approach 

The aim of this Section is to propose a redefinition of the semi-classical limit for gravity in 
minisuperspace which includes the (negative) energy of the gravitational fluctuations found in 
the previous Section into the HJ equation. This amounts to treat the gravitational fluctuations 
as an extra "matter" contribution, in much the same fashion as is usually done in perturbation 
theory around a fixed background in order to compute the back-reaction on the metric (see, e.g., 
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pO|). The main advantage of the BO decomposition ||, [| with respect to the latter approach is 
that we now derive such a description from a (supposedly more fundamental) unitary quantum 
theory (the WDW equation ( |2.1[ )) together with explicit conditions for the semi-classical approx- 
imation which should otherwise be deduced from external principles. In fact, this will give us 
(semi) classical trajectories Rf corresponding to the matter content ( Hm ) and the gravitational 



state / together with the consistency conditions discussed in Section 3.2 



In order to simplify the analysis from now on we shall consider one gravitational mode at 
a time and set ip = ip\ so that the energy of the gravitational fluctuations is given by E\ 
in Eq. ( 2.17 ). With the above restriction, Eq. ( |2.12j ) in the classical adiabatic approximation 
i?<i? becomes 



iKhip c P ^ = ip c 
oR 



-- 2 P>--Ri + R { H M 



2 A 2 



fx 



(3.1) 



The term of order (h/X) 2 survives in the semi-classical limit only provided one allows for very 
short wavelengths, such that ^ P /|A| does not vanish for Ti — > 0. This is just the analogue of what 
is required for the expectation value of the matter Hamiltonian (~ h/is), to wit £</, ~ fi. The 
condition |A| ~ £p, in turn, would refer to a fully quantum theory of gravity, if A is interpreted 
as a spatial wavelength, and one might prefer to place a ultra-violet cut-off A > i v for the values 
of A. We prefer to stick to a more euristic attitude and assign a physical meaning only to the 
energy E\, keeping it finite (and mostly small) throughout the computations. Of course one 
can always consider A ~ £ p as a (limiting) case of particular interest. Indeed, we will see in the 
next Section that one can obtain significant corrections induced by such modes in a way which 
is phenomenologically acceptable within the semi-classical treatment. 

One observes that the factorization of the wave-function ip\ (gravitational state) into f\ 
(fluctuations) and a specific ip c (classical part) is not forced by Eq. (|3.1|) or any other equation 
following from the WDW equation ( |2.1| ). In Refs. |^, |5], |9| it was rather determined by the 
implicit assumption that ( Hm ) is the dominant contribution in the semi-classical limit. This 
physical assumption takes mathematical form in the condition ( p. 22 ) which leads to the definition 
of the classical momentum P c in Eq. (|2.6|). However, since in the last Section we concluded that 
there are times at which E\ ~ ( Hm ) for every A, this is clearly contradictory and one should 
instead treat E\ as a source for the dynamics of R on the same footing as ( Hm ) ■ This can be 
achieved straightforwardly by introducing the modified momentum 



Px 



1 dR^ 



1 



2kR x {H m \ 



eR\ 



(3-2) 



A more formal way to derive this result is by defining a new WKB wave-function t/'a peaked 
on a modified trajectory R\, parameterized by a time variable rj\ {dr)\ = r] c dr] c ), such that 



l H 

e 



^P~x 



PxdRx -r 



(3-3) 



Then, upon substituting into Eq. fl2.2|) , in the classical adiabatic approximation R -C R and 
neglecting the RHS, one obtains 



-P 2 R 2 + R(H M ) o 

2 2/c X ' 2 A 2 



fx ■ 



(3.4) 



8 



The LHS vanishes identically, since —Kip\ P (df\/dR) = df\/dr]\ = 0. Therefore the RHS gives 
the modified HJ equation 



Ml 
p - 



(3.5) 



which is Eq. d3 



3.1 Examples 

The above Eq. (|3.5|) will now be solved in the quantum adiabatic approximation ( R~m ) = 
constant and for the three values taken by the parameter e for the purpose of showing explicit 
results. However, we emphasize that the latter approximation is not essential for the general 
formalism and is assumed just because it allows to carry on the computation analytically. 

Negative curvature 

For e = — 1 the velocity R\ vanishes at 



R 



± 







, 1± \ 





(3.6) 



thus one expects a turning point at R\ = R x (R x < is unphysical). The latter reduces to 
the turning point R~ = for ^/|A| — > 0. Upon setting R\(0) = R x , the modified trajectory is 
given by (rj = rj\) 



^ 1 + iv|^ cosh7? - 1 



(3.7) 



and the solution R c given in Eq. (^j|) is recovered as Roo in the limit ^/|A| — > (R x —> 0). In 
the opposite limit, |A|/£ P — > 0, R x diverges and the trajectory eventually reduces to a point. 



Flat space 

For e = the modified trajectory is given by 



f 2 Pj. NjJ 2 
„ > H_ , iv 0> 2 

2A2 21, V ' 



(3. 



with a turning point at R\(0). The solution R c in Eq. ( |2.9| ) is recovered in the limit ^ P /|A| — ► 
(R\(0) — ► 0) as for e = — 1. Also, the opposite limit behaves the same as for negative curvature. 



Positive curvature 

For e = +1 there are two turning points at 

IJ2 



i ± 



\ 



N 2 A 2 



(3.9) 
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provided the square root is real, that is |A| > 

Rt 



R 



A 



R- 



and the maximum Rt 



the modified solution is 



2N A £ 2 p 



As before, for ^/|A| — > the minimum 
(f,tp/£^. Again, upon setting -Ra(O) = R 



A > 



\ 



1 



2 Yl cosr ? 



iVl A 2 



(3.10) 



Thus, the effect of the extra term in the HJ equation is to make R\ oscillate between a minimum 
value R^ which is shifted above zero and a maximum value which is below the turning point 
R£ . The shifts vanish and the solution R c given in Eq. (2.9) is recovered in the limit ^/|A| — ► 0. 
At the opposite limit stands the case |A| = i^/N^ for which the amplitude of the oscillation 



Rt - r: 



0. 



3.2 Consistency of the approximations 

The trajectories R\ displayed above are consistent only when three different approximations 
hold simultaneously: 

1. 



The RHS's of Eqs. ( j2.2|) and ( |2,4D must be negligible. This is the only relevant approxima- 
tion which must hold for the general formalism as developed in this Section to apply. Since 
the RHS of Eq. (|2.2j) amounts to the expectation value of a gravitational operator on the 



matter state x, it is not substantially modified by the new definition of in Eq. (2.11) with 
respect to the state ( |2.5|) (one only expects corrections because x W1U evolve differently 
with the new scale factor R\). Hence, we refer the reader to Refs. [jll], [l2|, 13, 15, 16|, 
where this approximation has been studied at best for the case in Eq. (^5|). 

However, in order to ensure that the corrections we have computed are significant, we now 
need to check that the RHS of Eq. ( [2.2| ) is negligible with respect to the new term \E\\. 



Taking the estimate for the RHS as given in Appendix B of Ref. [11 



RHS 



^4 N 2 
R 3 J V 



we obtain the condition 



(3.11) 



RHS < \E X \ 



R^N S \ 



(3.12) 



where we keep on singling out the particular value A ' 
given in the previous Section and find that, for e 
satisfies Eq. ( |3.12 ) provided 



i p . We can now take the trajectories 
0,-1, the turning point Rx~e p (fy 



(3.13) 



while for e = +1 such a condition would be unphysical and is never met. This means that, 
for e = +1, the semi-classical approximation breaks down before any rebound occurs and 
quantum fluctuations prevail at small R so that the superposition among different semi- 
classical trajectories cannot be avoided. Strictly speaking, the latter conclusion extends to 
all values of e, since the WKB approximation breaks down near the classical turning points 
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This actually turns out to be a blessing in disguise, since, in the case of the collapse, the 
rebound of the sphere at a finite radius inside the Schwarzschild radius would eventually 
violate causality |29[| . 

A more complete analysis of the dynamics when the RHS's are not negligible is currently 
under study f2£|. It is expected that one can still determine the evolution of the system 
in this non-classical range (at least numerically), although the geometrical interpretation 
is then lost (see, e.g., Ref. 



The quantum adiabatic approximation ( Hm ) constant. As we mentioned previously, 
this approximation is not essential for the development of the general formalism and was 
actually relaxed in Refs. [13, 15, [l(|. However, it is only for ( Hm ) constant that solutions 
can be computed analytically. Again we refer to Refs. [11, 12, 13, 15, 16] for a detailed 



analysis. For the case of the gravitational collapse, one obtains the condition [11] 



r - r(0) > 1$ , 



(3.14) 



where r = p s R\(r]) is the areal radius of the sphere and r(0) is the turning point p s R\(0). 

The classical adiabatic approximation R\ <C R\. Despite the terminology, this approxi- 
mation is not related to the previous one, although, to some extent, it is inessential for the 
general formalism developed here. However, it is necessary for Eq. ( [2.14 ) for / to simplify 
to the tractable form ( |2.15 ) and leads to a new condition which turns out to be not really 
restrictive for the example worked out here. In fact, it is easy to see that for the proposed 
solutions the ratio 



Rx 
Rx 







(3.15) 



for A -> (e 



-1,0) and A 



(the minimum allowed value for e = +1). 



From the above considerations one thus concludes that the break-down of the semi-classical 
approximation might occur at relatively large values ~ R\(0). In the following Section, such 
turning/breaking points in the gravitational collapse will be regarded as significantly different 
from the point-like singularity only if they are bigger than the Compton wave-length of the 
particles. 



4 Applications 

As described in Section |2], the contribution of gravitational fluctuations incorporated in the 
theory is of quantum origin. It is generally taken for granted that gravity in the world we 
can test is classical, which leads one to assume the energy stored in quantum gravitational 
fluctuations is negligible at the time 770 when measurements take place. To be more precise, let 
us introduce the ratio a 2 (r/) = \E\\/(Hm) between quantum gravitational fluctuation energy 
and matter energy and assume ao = 01(770) <C 1. From Eq. (|2.17|) this definition can be used to 
express A as 

\ 2 - — — — E-t- (A 1 \ 
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where R = R\(r]o)- 

The second important issue is whether the inclusion of E\ in the HJ equation leads to 
observable effects, that is, one will have to check when (if ever) deviations from the standard 
trajectories R c are physically significant in magnitude. 

In order to clarify the above points, we now specialize the very simple analytic solutions 
found so far to two models. In so doing we do not expect our results to be definite answers to 
any basic physical question in either cases, however, we believe they give hints as to the possible 
relevance of the predicted effects. 



4.1 Cosmology 

In the cosmological case one takes 770 equal to the (conformal) age of the Universe, so that the 
energy stored in the gravitational fluctuations is totally negligible today. However, this does 
not prevent a 2 ~ 1/-Ra to be comparable with one or bigger at very early stages. The key 
observation is precisely that the present scale factor of the universe, Ro, is related to the initial 
(minimum) scale factor R(0) by 

al = a\0)^. (4.2) 

One also recalls that in the RW metric the spatial distance between two points arbitrarily set 
at p = and at p = p^ is given by (pd < 1 for e = 0) 

. = */"-*. (4.3) 



Open Universe 

When the spatial curvature e = — 1, from Eqs. ( |3.6| ) and ( ^4.1[) one obtains 



R(0) = R- = N^ 
On using Eq. fl4.2|) one finds 



10 qN+ 



(4.4) 



R(0)=2N 4 ,f (a 2 (0)-l) , 



(4.5) 



so that -R(O) > implies a 2 (0) > 1 and quantum gravitational fluctuations must dominate the 
early stages in order to have a start at non-zero scale factor (this is due to the gravitational 
potential contributing with the same sign as ( Hm ) in the HJ equation for e = —1). 

It is now interesting to relate «o to physically meaningful quantities. For instance, the 
relative difference As = Sf — s c between spatial distances measured when R = Rx^t/q) and 
R = R c (r]o) is 



Asq 

so 



1 + 2q 2 



Rq 



(4.6) 
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where we used Eqs. fl2.9| ) and (3/7) in the limit coshr/p 3> 1 (Asp > since R\ ~ i? c + i? A ). 
Hence Eq. ( [4.4| ) becomes 



i2(0) ~ iV> / 



A.s 
so 



(4.7) 



If one takes for Asp/sp the accuracy with which distances are measured in the present Universe, 
the above equation gives us a (very rough indeed) estimate of the maximum value for the initial 
scale factor which cannot be ruled out by present measurements. 

A further consequence of having R(0) > is that two points in space were causally discon- 
nected at r\ = if their distance s(0) 3> 1$, that is, from Eq. 



Pd > Pc = sinh 



•so 



(4.8) 



In the flat case e = 0, a 2 (0) = 1, as follows from a quick inspection of Eq. (3.5), and the 
expression of the initial scale factor simplifies to 



R(0) = al R 



Asp 
«o 



Rq 



(4.9) 



Then the causal comoving radius p c = ■£ 



so 



Closed Universe 



For e = +1 we have shown in Section 3.2 that our corrections are less relevant, however, we also 
consider this case for completeness. On taking R ~ R^, one has 



Ro^2N 4 



1 - a; 



o > 



(4.10) 



which holds for < 1 (the gravitational potential in the HJ equation is now opposite in sign 
to (Hm ))• The relative difference in distances is 



Asp 

«o 



-a, 



! 



(4.11) 



where we used Eqs. ( |2.9D and ( |3,10| ) in the limit cosr/o ~ — 1 (Asp < because the maximum of 
R\ is shifted down with respect to the maximum of R c ). Putting the pieces together gives 



and s(0) > 1$ if 



R(0) 



P^> Pc = sin 



As 



so 



Rn 



Ro 



so 



Asp 



Since p < 1, the latter condition can be satisfied only if 



Asp 

so 



(4.12) 



(4.13) 



(4.14) 
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No need to say the present model is too simplified to take the expression for the comoving 
radius p c seriously as a prediction for the scale of Cosmic Background Fluctuations or related 
cosmological quantities. In fact, there is no inflationary stage and regions outside of p c will 
eventually come into the causal cone after a finite (short) time due to the slow expansion of 
the scale factor in the FRW models. The situation might change in case one considers a more 
realistic description. Further, we recall that the existence of a minimum scale factor -R(O) > 
is a basic ingredient of Pre-Big-Bang Cosmology (see, e.g. and Refs. therein), in which case 
(f) should be identified with the homogeneous mode of the dilaton predicted by the low energy 
limit of string theory [32]. 



4.2 Gravitational collapse 

For the case of the collapsing sphere of dust the above framework is inverted since now \E\\ 
increases along the classical trajectory. Thus, although one starts with «o < 1 so that the energy 
of the quantum gravitational fluctuations is totally negligible, when the singularity is approached 
the gravitational fluctuations induce an effective quantum pressure which slows down the collapse 
and causes the break-down of the semi-classical approximation at a finite radius of the sphere 
r(0). A very important observation, already mentioned at the end of Section |3|, is that r(0) is 
physically distinguished from the singularity p s R c = only if it is bigger than the Compton 
wavelength 1$. 

Before proceeding, it is useful to recall that the Schwarzschild radius of the sphere rjj = 2M 
where the ADM mass parameter is 

i 2 

M = p 3 s N <i ,f, (4.15) 

regardless of the value of e. We will assume M/k is the mass that is measured for astronomical 
objects, although it generally differs from the proper mass (h/i^) J Ps p 2 dpj \J\ — e p 2 . Fur- 
ther, p s is related to the geodesic energy parameter £ of the trajectory r s = p s R of the radius 
of the sphere in the outer Schwarzschild space-time with mass parameter M by £ 2 = 1 — e p 2 
(— 1 < £ < 1 for bound orbits, £ > 1 for unbound orbits) [26[ |. 



Scattering orbits 

For £ > 1 one can choose the starting radius of the sphere ro = p s Ro is any value greater than 
r#. Then the sphere will bounce in correspondence to i?^ at 



r(0) 



M 



l + 2{£ 2 -l)a 2r f i 



1 



£ 2 -l 

As mentioned above, r(0) must be greater than 1$ to be physical, that is 

[r H + (£ 2 -l)U 



For £ 2 — 1 small one can expand the square root in Eq. ( 4.16j ) and obtains 

r(0) ~ «q r o , 



(4.16) 



(4.17) 



(4.18) 



14 



with the condition a 2 , 3> I^/tq. The above result is exact for £ = 1, in which case p s is arbitrary 
and r = (a 2 , vq) + p s R c . In the opposite limit, £ S> 1, one obtains 



r(0) ~ ^ ^r r H 



(4.19) 



with a 2 , ^£ 2 l\jr H r . 

In all cases (e = —1,0), r(0) is bigger than r# provided Oq > r#- £ 2 /r , or, using Eq. (4.1) 
and ((LTD, the comoving radius is given by 



2M 



(4.20) 



for the limiting case A 



Bound orbits 

For £ 2 < 1, on setting ro = p s i?^, the choice (4.1) and the trajectory ( |3.10 ) give 



ro 



l-£ 2 



(l-at) . 



Then, the radius at which the sphere bounces is given by 

M 



r(0) 



l-£ 2 



1_ Jl-Aa 2 (1-a 2 ) 



(4.21) 



(4.22) 



and one concludes that Eq. ( 4.1£| ) holds for — 1 < £ < 1, in which cases r(0) 3> provided 



<C 1 and r(0) > rn for r^/ro < -C 1. For the limiting value A ~ ^ p this means 



9 2M 



(4.23) 



which lies outside the allowed range of p s and the breaking point cannot thus be bigger than m 
for e = +1. 

It is worth noting that the minima obtained above are fairly generic in that they do not 
depend on the detailed structure of the sphere nor on the specific form of the quantum grav- 
itational fluctuations (i^, t v and A do not appear explicitly). Furthermore, the turning point 
at which the semi-classical approximation breaks can be rather big and, from Eq. ( 4.1S| ), one 
cannot exclude it occurs at a radius comparable with the error with which one measures ro. Of 
course, to model a star with a sphere of dust ignores (among the rest) the crucial role played by 
the pressure in keeping the star in equilibrium and contrasting the collapse itself. Hence, it is 
clear that the actual value of r(0) could be significantly different from the one estimated here. 

We conclude by mentioning an independent argument which supports the possibility of 
having the kind of pressure emerging from the quantum fluctuations discussed in this Section. 
In fact, besides the WD W equation, one has the conservation of the total energy of the system, 



namely its ADM mass [33 



M = p% UH M )-E f )~pi(H M (r l0 



(4.24) 
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Since ( Hm ) increases in time, due to non-adiabatic production of matter particles [jl^, [l5| (an 
effect totally ignored in the present notes), Eq. ([4.24 ) requires that either p s decreases or Ef 
increases (or both effects take place). The first case amounts to quantum jumps to classical 
trajectories with geodesic energy closer and closer to £ = 1 which would act as a semi- classical 
attractor ||34|| . The second case would imply that, although one can start with a state in which 
only matter modes are present (as appears sensible for a sphere of large initial radius), the price 
to pay for preserving the classical dynamics is the generation of gravitational perturbations in 
an amount such that the total energy (along with £ ) is conserved. We remark that the spherical 
symmetry assumed for the model would prevent these perturbations from propagating in the 
external vacuum as gravitational waves (Birkhoff 's theorem). 



5 Conclusions 

We have generalized the BO approach to the WDW equation in FRW minisuperspace || in 
order to include homogeneous quantum gravitational fluctuations around the WKB trajectory. 
In a standard approach, such fluctuations are treated as perturbations of a fixed background 
and satisfy a Schrodinger equation, whose solutions in (double) adiabatic approximation were 
displayed both in Lorentzian and Euclidean space. One then realizes that these solutions sig- 
nal a possible instability, since their energy can grow without bound for small R. The latter 
result suggested that the semi-classical limit had been incorrectly identified. Then a second 
approach was proposed in which the semi-classical limit includes the back-reaction of quantum 
gravitational fluctuations on the metric from the start, in much the same fashion as in quantum 
field theory on curved back-ground one replaces the classical energy-momentum tensor of matter 
with the expectation value of the corresponding quantum operator (for matter and gravitational 
waves). This in fact led to different classical trajectories with a possible non- vanishing minimum 
size for a FRW Universe or a collapsing body at which the semi-classical approximation breaks 
down and the metric does not admit a description in terms of classical variables. 

Although the FRW minisuperspace model is too simple to be realistic, one learns that the 
role played by quantum gravitational fluctuations, at least, should not be overlooked when 
considering self-gravitating matter. The guiding analogy is the treatment of infrared divergences 
in quantum field theory. In that context one encounters diverging quantities when studying, 
e.g., the Bremsstrahlung from an electron moving in the external electric field of a nucleus. 
At the classical level (the tree level of the quantum theory) the transition amplitude diverges 
for vanishing energy of the emitted photon. This problem is cured by adding the (diverging) 
one-loop contribution and noting that experimental measurements would not distinguish the 
final state of the electron with energy E from any state with energy E — AE if AE is smaller 
than the precision A of the apparatus. Then one has to sum over all the (tree level) emissions of 
(soft) photons of energy AE < A and obtains the counter-term which precisely cancels against 
the one-loop diverging term. Perhaps one can rephrase the results obtained in these notes by 
saying that the inclusion of (soft) quantum gravitational fluctuations with energy smaller than 
the precision with which we measure the energy of matter seems to cure the singularity in 
the density distribution of matter which develops at the classical level according to General 
Relativity. 
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